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Abstract 

We present the results of a detailed study of energy correlations at steady state for a 1-D 
model of coupled energy and matter transport. Our aim is to discover — via theoretical argu- 
ments, conjectures, and numerical simulations — how spatial covariances scale with system 
size, their relations to local thermodynamic quantities, and the randomizing effects of heat 
baths. Among our findings are that short-range covariances respond quadratically to local 
temperature gradients, and long-range covariances decay linearly with macroscopic distance. 
These findings are consistent with exact results for the simple exclusion and KMP models. 
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Introduction 

Transport processes, such as heat flow through a conducting medium in contact with unequal 
heat reservoirs, are intrinsically nonequilibrium phenomena because of the presence of nonzero 
currents J51[T7l- A central problem in nonequilibrium statistical physics is to explain how such 
large-scale, macroscopic processes arise from complex microscopic interactions. The most basic 
questions are perhaps those of mean profiles of quantities of physical interest and their responses 
to external forces. Nonequilibrium steady states, on the other hand, are well known to be char- 
acterized by large fluctuations, and among the simplest measures of fluctuations are temporal and 
spatial correlations. 

In this paper, we present a systematic study of spatial correlations at steady state for a class 
of 1-D stochastic models called Random Halves Models. A detailed description of this class of 
models is given in Sect. 11.11 Very briefly, a linear chain of open cells is connected to two unequal 
heat baths, which inject into the chain tracer particles at characteristic rates and energies. Energy 
storage devices are systematically placed throughout the chain to mark local energy levels. Their 
contents are redistributed by the tracer particles as they move through the chain. These models 
were introduced in |[9l as stochastic idealizations of certain mechanical models lfT6l[T8Tl . 

Our study is based on a combination of analytic arguments and numerical simulations. We 
believe random halves models are excellent candidates for this method of investigation for two 
reasons. First, their dynamics are richer and more complex than rigorously-understood models 
such as the simple exclusion [6] and KMP models lfT3l . as random halves models have two trans- 
ported quantities (energy and matter) and highly nonlinear interactions involving two different time 
scales. These features make a purely analytical study more difficult. Second, all of the forces acting 
on this system are clearly identified. This is seldom the case in more realistic physical models. 
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Our main results can be summarized as follows. Given boundary conditions, i.e., the tem- 
perature and injection rate of each bath, we let Si denote the stored energy at the ith site, and 
Cov N (Si, Sj) the covariance of Si and Sj in a chain of length N. Our first finding is Cov^ (Si, Sj) ~ 
jr for i ^ j, which leads us to consider the following two functions describing the covariances of 
stored energies at microscopic and macroscopic distances: 

C(x) = lim TV • Cov N (S[ xN ],S [xN ] +1 ) , x £ (0, 1); 
C 2 (x,y) = lim N ■ Cov N (S [xN] , S [yN] ) , x,y e (0,1), x^y. 

We show that C(x) has the form 

C(x) « bc (x) ■ A(k(x)) ■ T\xf (1) 

where the right side is to be interpreted as follows: (p hc (x) is a measure of effective distance to 
the baths, a rough approximation of which is (p bc (x) ~ 4r(l — x); T(x) is local temperature, and 
T'(x) is temperature gradient; k(x) is local particle density, and A(k) — > as k — > oo. Note the 
quadratic dependence on the local temperature gradient. For long-range covariances, we find that 
for fixed x, the function y i— > C 2 (x,y) is continuous but not differentiable at y = x, decreasing 
roughly linearly as \y — x\ increases. 

Our numerical results are consistent with exact analytic results on the simple exclusion and 
KMP models (6l|71|20l[3]|. We extend the existing picture to a setting of coupled energy and matter 
transport, with features not present in these two previously studied models. For example, we 
show that energy covariances at microscopic distances respond quadratically to local temperature 
gradients (which vary along the chain; see [0). There is a second transported quantity, namely 
matter in the form of particles, and energy covariances are shown to be inversely related to local 
particle density. Throughout the paper, we provide cross-checking numerical evidence for the 
phenomena we identify, and venture to give a physical interpretation whenever we can. See dlU 
[TT1 [14l H5J for other relevant works on this topic. 

While most of the results presented here are specific to 1-D, they also serve as a basis for direct 
comparisons with higher dimensions. A corresponding analysis of random halves models in 2 and 
3-D is underway. The results will be reported in a separate paper. 

Note on simulations. All of our numerical results are obtained via direct simulation. That is, 
our computer programs faithfully implement the dynamics of the random halves model described 
in Sect. |1.1[ and expectation values with respect to the invariant measures are computed as time 
averages over long trajectories. Some relevant numerical issues are discussed in Appendix A. 

1 The Random Halves Model 

1.1 Model Description and Physical Interpretation 

A class of models of nonequilibrium phenomena is introduced in [I9J. In each model, there is a ho- 
mogeneous conducting medium represented by a linear chain of N identical "cells" with stochastic 
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Figure 1: The rotating disc model. The "random halves" model studied in this paper is a stochastic version 
of this model. 

heat baths coupled to the ends of the chain. Within each cell, there is a stored energy which charac- 
terizes the "temperature" at that location. Matter (in the form of tracer particles) and energy (in the 
form of tracer kinetic energy) are injected into the system by the heat baths; they are eventually ab- 
sorbed by the heat baths. Tracer particles interact with the local system at each site, redistributing 
the stored energies as they move about in the chain; they do not interact directly with each other. 
This general framework was introduced as an abstraction of the mechanical models in lfT6l[T8l . It 
encompasses both "Hamiltonian" models with conservative deterministic dynamics (the heat baths 
being the only sources of randomness), and stochastic models, in which the microscopic dynamics 
are defined by conservative stochastic rules. 

For concreteness, we begin by describing a slightly simplified version of the mechanical models 
studied in [fT6l[T8l . In this rotating disc model, each cell contains a disc nailed down at its center, 
about which it rotates freely (see Fig. [T]). Whenever a tracer collides with a disc, it exchanges 
kinetic energy with the disc via a deterministic rule, e.g. the angular momentum of the disc may 
be interchanged with the tangential component of the tracer's momentum. When a tracer collides 
with a cell wall, it reflects elastically. Here, the stored energy at each site is the rotational energy 
of the disc. 

The random halves model studied in the present paper are stochastic idealizations of mechani- 
cal models like the rotating disc model described above. Nonequilibrium energy and particle den- 
sity profiles of random halves models have been analyzed in [|9]|. The rest of this section reviews 
relevant parts of that paper. 



Precise description of the Random Halves Model. There are iV sites, labeled 1, 2, ■ • • N, with 
tracers moving through the chain. At each site there is an abstract energy storage tank. We let Si(t) 
denote the amount of energy in the tank at site i at time t. The microscopic dynamics are defined by 
the following rules. Fix 5 > 0. Each tracer is equipped with two independent exponential clocks. 
Clock 1, which signals the times of energy exchanges with tanks, rings at rate \Je(t) /8, where e is 
the kinetic energy of the tracer. Clock 2, which signals the times of site-to-site movements, rings 
at rate \/e{t). Suppose a tracer is at site % and one of its clocks rings. Then instantaneously: 

(i) If Clock 1 rings, the tracer energy e and the stored energy S\ are pooled together and split 
randomly. That is, the tracer gets U ■ (e+Si) units of energy and the tank gets (1 — U) • (e+Si), 
where U E [0, 1] is uniformly distributed and independent of all other random variables. 
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(ii) If Clock 2 rings, the tracer leaves site i. It jumps with equal probability to sites i ± 1. A 
tracer entering site or site N + 1 exits the system forever. 

All tracers originate from and eventually exit to one of the heat baths. Each heat bath injects tracers 
with independent, exponentially-distributed energies into the system. The left bath injects tracers 
with mean energy Tl into site 1 at an exponential rate of pl, and the right bath injects tracers with 
mean energy Tr into site N at a rate of pR. Tracers in the system are indistinguishable. 

Physical interpretation and remarks. It is natural to think of the parameters Tl and Tr as 
temperatures. This leads us to define the temperature T at site % in an iV-chain to be T := E(5j), 
where the expectation is taken with respect to the invariant measure of the iV-chain. The injection 
rates p L and p R can be rewritten in terms of chemical potentials; we do not pursue this analogy 
further. 

The Hamiltonian model contains a small length parameter which does not appear in the stochas- 
tic model, namely a length t which measures the size of the cell. In the stochastic model, we set 
t = 1. It is useful to keep this in mind in dimensional analysis. 

There are two time scales in our system, one associated with the local dynamics at each site 
and the other the movement of tracers along the chain. The ratio of these two time scales is S. For 
example, 5 < 1 means that on average, tracer-tank energy exchanges occur much more frequently 
than site-to-site movements of tracers. Steady-state macroscopic profiles such as temperature and 
tracer density do not depend on 5, but 5 can have a significant impact on the numerical values of 
spatial and temporal correlations. In this paper, we have chosen to simplify matters by fixing 5, 
which throughout the paper is set equal to i. 

At any moment in time, the number of tracers at each site can vary from to oo. Observe 
also that the interaction in this simple model is highly nonlinear, even though there are no direct 
tracer-tracer or tank-tank interactions: tracers at the same site exchange energy via interacting with 
the tank, and all actions — including energy exchanges and site-to-site jumps — take place at rates 
proportional to the "speeds" (the square roots of the kinetic energies) of the tracers at that moment 
in time. This is essential if these stochastic models are to mimic the behaviors of their Hamiltonian 
counterparts. 

1.2 Invariant Measures at Equilibrium 

At equilibrium, i.e., when the left and right baths have equal temperatures and injection rates, 
the invariant measure of the random halves model is known explicitly. To give its density, we 
need a little bit of notation: the state of a single cell in which there are exactly k tracers is spec- 
ified by (S, {ei, e 2 , • • • , e^}) where S > is the stored energy, i.e., the energy of the tank, and 
{ei, e^i ' • • j e k} is an unordered set of k positive numbers representing the k tracer energies^. Let 
denote the set of all possible states of a single cell with exactly k tracers present. The state 
space for a single cell is then the disjoint union f2 = Ufcf2 fc . 

'it is an unordered set because the tracers are indistinguishable. 
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Proposition 1.1. /[9|/ Let T L = T R = T, and Pl = Pr = P- Then the unique invariant probability 
measure of the N -chain is the N-fold product 

p N = p T ' p x • • • x p T ' p 

where p T ' p is the measure giving the statistics within each cell. The measure p T ' p is defined by: 

(i) The number of tracers present is a Poisson random variable with mean k = 2p^Jn/T, i.e., 

k 

p T > p (n k ) = p k := ^-e- K , k = 0, 1, 2, ... . (2) 
k\ 

(ii) The conditional density of p T)p on Vt k is c k a k ({ei, ■ ■ ■ , e k }, S) where 

a k ({e u . . . , e fc }, S) = 1 e -^+- + e k+ s) (3) 

V e l • • • • " e k 

k /2 

here j3 = 1/T, and c k — j3k\(j3/ir) is the normalizing constant. 

Proposition 11.11 shows that at equilibrium, the stored energies at different sites are entirely 
uncorrelated. A generalization of this result is given in Proposition 5.1. 

1.3 Macroscopic Equations of Nonequilibrium Steady States 

We now fix arbitrary boundary conditions T L , T R , p L , and p R , and give the equations of temperature 
and tracer profiles in nonequilibrium steady states. 

Assumption 1. We assume that for each N, there is a unique invariant probability measure p^ to 
which all initial data converge. 

Because the energies are unbounded, a tightness argument is needed to guarantee existence. 
Uniqueness should be straightforward. In what follows, the word "mean" refers to expectation 
with respect to pn- We identify each site i with the point X; = j^?— in the unit interval, and think 
of x = and x N+ i = 1 as the locations of the baths. For i = 1, • • • , N, define 

p(xi) = mean number of jumps from site i to site i + 1 per unit time 
= mean number of jumps from site i to site i — 1 per unit time; 

q(xi) = mean energy flow from site i to site i + 1 per unit time 
= mean energy flow from site i to site i — 1 per unit time. 

That is to say, the mean number of jumps out of site i per unit time is 2p(xj); half go to the right 
and half to the left, and so on. The quantities p and q are well-defined in steady state and have very 
simple behavior: 

Lemma 1.1. /[9]/ With p and q defined as above, we have, for x = x iy i = 1, • • • , N, 
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(i) p(x) 



= Pl + (pr - Pl)x; 



(ii) q(x) = p L T L + (prT r - p L T L )x. 

As N — > oo, these functions converge (trivially) to linear or, more accurately, affine func- 
tions on (0, 1). Note that p' and q' are the steady-state tracer and energy currents, respectively. 
These currents by themselves, however, do not determine steady state profiles such as those for 
temperature and tracer density. The following assumption is used to deduce such information: 

Assumption 2. (A version of LTE) For 1 < i < N, let p Nji denote the marginal of p^ at the 
site i. We assume that for every x G (0, 1), Pn,[ x n] converges as N — > oo to p T,p ^ for some 
T = T(x) > 0. A tightness condition for all the Pn,[xN] is also assumed. 

Let ki denote the number of tracers present at site i. Recall that Si is the stored energy. 

Theorem 1. /[9]/ Let arbitrary boundary conditions T L , T R , p L , p R be given. Under Assumptions 1 
and 2, 

T(x) = lim E m (S [xN] ) and k(x) = Jim E m (k [xN] ) (4) 
are well-defined and are given by 

T(x) = g} an, «(,) = 2^-^. (5) 

See |H for details. 



2 Global Fluctuations and Pair Correlations 

This paper concerns the spatial energy covariances at steady state of the random halves model 
described in Sect. ll.lt More precisely, let boundary conditions T L , T R , p L and p R be specified, and 
consider an iV-chain. The quantities of interest are 

Cov^S,, S,) = E^SiSj) - E m (Si) E m (Sj), l<i,j<N, j. 

We begin by examining how these quantities scale with N. 

2.1 Scaling of Total-Energy Variance 

For an iV-chain with boundary conditions Tl, T r , pl and p R , we consider the total energy variance 

N 

V(N; T L , T R , p L , p R ) = Var ( ^ S?) . (6) 

i=l 
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Figure 2: Variance of total energy as function of system size. The upper curve is V(N), the lower curve is 
V Q (N). The boundary conditions are T L = 10, T R = 100, p L = 20, and p R = 10. 

Keeping the boundary conditions fixed, we treat V = V(N) as a function of N, and observe that 
V = Vq + Vi where 

N 

K, = ^Var(5 i ) and ^ = £ £ Cov^, S,). (7) 

i=l i jy^i 

By Theorem 1 , 

l r 1 

lim -V (N) = / T(x) 2 dx. (8) 

That is, V (N) ~ -BoA" for iV > 1 with fio = /„* T(x) 2 rfx. 

When T L = T R , V x = (see Sect. so that V(JV) = Vo(A^) + Vi(N) ~ AT. We expect 
this to be true when 7^ Tr on physical grounds, which leads to the question: is it true that 
Vi(N) ~ B\N, and if so, what is the sign of B\l Two sets of numerical simulations are performed 
to resolve this. In one, we compute V\ — V — Vq directly as a function of N . In the other, we start 
from equilibrium, and hope to observe a quadratic response as a temperature gradient is introduced; 
the sign of the quadratic term is then the sign of B\. The reason we expect a quadratic (as opposed 
to linear) response is symmetry: the temperature gradients T R — T L = AT and T R — T L = —AT 
clearly lead to the same energy variances and co variances. 

Simulation 1. Fix bath temperatures Tl,T r and injection rates Pl,Pr and compute V(N) and 
Vq(N) for increasing N. 

The results, plotted in Fig. [2l show that 

V(N) ~ BN with B> B >0, (9) 
i.e. Vi(N) ~ B X N for some B t >0. 

Simulation 2. Fix N, pi = p R , and a number m > 0. Compute V and Vq for various pairs 
(Tl, T R ) chosen so that TL ^ Tfl = m, and investigate V\ — V — Vo as a function of AT = T R — Tl. 
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The results show that V\ depends quadratically on AT with a strictly positive coefficient. Taken 
together, these two simulations suggest that for p L = p R and fixed mean temperature, 

V 1 (N; T L , T R , p L , p R ) cc (T R - T L f ■ N , (10) 

at least when \T R — Tl\ is not too large compared to the mean temperature in the chain. We will see 
later that both the linear scaling of V\ with system size and its quadratic response to temperature 
gradient are consistent with the scaling of Cov N (Si, Sj). 

Note that one advantage of probing spatial correlations via global quantities such as V\ is that 
it can be computed more reliably than small, local quantities such as Covjv (Si, Sj). 

2.2 Pair Covariances 

We begin now to investigate the individual terms in the sum Vi = Ylj^i Cov A r(S' i , Sj). 

The two most basic characteristics of Cov N (Si, Sj) are its sign and order of magnitude. The 
fact that V\(N) ~ N suggests that many of the covariances are positive. Indeed, all of our nu- 
merical evidence points to Cov N (Si, Sj) > 0, as does the theoretical reasoning in the sections to 
follow. There is reason to remain cautious, however: Since Cov N (Si, Sj) = when the system is 
in equilibrium (Proposition 1.1), and can be zero even when there is a tracer flux (Proposition 5.1), 
one cannot conclude definitively via numerics alone that there are no strictly negative correlations. 
See the remark at the end of this subsection on why the nonnegativity of covariances in this model 
may be a delicate question. 

We proceed to an analysis of the order of magnitude of Cov N (Si, Sj), assuming in the heuristic 
discussion below that Cov n(S%, Sj) > for all i, j. 

We propose to decompose V\ into 

Vi = Vj.j where V lti = ^ Cav N (S i7 Sj) 
and reason as follows: 

(i) For i away from the two ends of the chain, V^ ~ 1 since these terms sum to V\ ~ and 
it is unlikely that they scale differently with N. (For i close to the boundary, V\ :i may be 
smaller due to the entrance of tracers with i.i.d. energies.) 

(ii) For each fixed i, the function j \—> Cov]y(Si, Sj) should decrease monotonically as \i — j\ 
increases: sites farther apart are expected to be less correlated because they "communicate" 
via longer and noisier "channels." 

These considerations imply that Covjy(Si, Si+i), with i away from and N, are among the larger 
of the ~ iV 2 terms in V\. Since Covjv(5'i, Si+i) is one of N terms in V\^, (i) above implies 

Cov N (Si,S i+1 ) > — . 

We point out that the order of magnitude of Cov N (Si, S i+ i) contains a fair amount of in- 
formation about the shape of the function j \— > Cov N (Si, Sj): if Cov N (Si, S i+ i) ~ j?, then 
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Figure 3: Scaling of Covtv with N. This plot shows N • Covn{S[n/2], Sjjv/2]+i) as a function of N. 
Boundary conditions are Tl = 10, Tr = 190, and pi = PR = 20 for the upper curve and pl = PR = 40 
for the lower one. 

Cov N (Si, Sj) ~ jj for a definite fraction of j (because these terms have to add up to ~ 1). This 
points to an extremely slow decay of Cov n(Si, Sj) with increasing distance between sites i and j. 
If, on the other hand, Cov N (Si, S i+ i) 3> jj, then the function j i-> Covn(Si, Sj), j ^ i, would 
have a very sharp peak at j — i ± 1. The question is: which scenario is the case here? 
To resolve this issue, we perform the following numerical simulation: 

Simulation 3. For fixed boundary conditions, compute Covn(S[n/2], S[n/2]+i) for a range of N. 

Two sets of results are plotted in Fig. [3] They (and other data sets not shown here) show that 

C0V N (S[ N / 2 ], S[ N / 2 ]+i) ~ — • (11) 

As explained above, one may infer from this that all pair covariances Covjy(Si, Sj) = 0(1/ N), a 
fact confirmed by many subsequent simulations (see e.g. Fig. [T2)) . 

Remarks on sign of Cov N (Si, Sj). As far as we know, existing analytic techniques for proving 
or disproving positivity of correlations do not apply to the random halves model. This, in part, 
is due to the fact that in our model, mechanisms conducive to both positive and negative corre- 
lations are at work, and the actual covariance reflects a balance between these two tendencies. 
For instance, a large, upward surge in the energies of incoming tracers has the effect of raising 
some of the stored energies above their expected values. It is believed that such excursions lead 
to long-wavelength fluctuations with slow relaxation times, resulting in positive correlations lUOl . 
However, upon interacting with one of the tanks, a tracer that acquires a higher- than-normal frac- 
tion of the energy is likely to move more quickly to a neighboring site and to interact with the 
tank there, possibly transferring its energy to the tank in its new location. Such a phenomenon cre- 
ates negative correlations between neighboring sites. In equilibrium, these two tendencies balance 
perfectly, giving zero covariance (see Proposition ll.il ). Further numerical evidence in support of 
nonnegative covariances is given in Sect. 15.11 
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3 Covariance Structures "at N = oo" 



The purpose of this short section is to make explicit the objects to be studied in the rest of this 
paper and the assumptions under which we plan to operate. 

The discussion in the Sect. 12.21 points to the following functions: Given boundary conditions 

Tl, T r , p L and p R , and x, y G (0, 1), x ^ y, define C(x) and C 2 (x, y) by 



if these limits exist. 

Supposing these functions are well defined, it is natural to ask (i) how C(x) depends on local 
thermodynamic quantities such as T(x), k(x) and their gradients, and (ii) if these local quantities 
completely determine C(x). We present below an answer to the second question; it is one of a 
number of possible answers. Question (i) is addressed in the next two sections. 

Given x G (0,1) and (T L ,T R , p L , p R ), Theorem 1 gives a quadruple (T*, p*, Tl, p[) where 
T = T(x), p* = p(x), Tl = T'(x) and p'^ = p'. Conversely, given x G (0, 1) and a quadruple 
(T*, p*,Tl, p'^), it is easy to check that one can solve the equations in Theorem 1 backwards and 
find the corresponding boundary conditions. This one-to-one correspondence between boundary 
and local conditions defines a function T with 



Since we work only with finite chains, it is necessary to have a version of T before the infinite- 
volume limit. Given (T L , T R , p L , p R ), x and N, we let % = [xN], and let T, p h AT and Api be 
as follows: T is the mean temperature at site i, pi is the "injection rate" (see Sect. II .31 for precise 
definition), ATj = T i+1 — T and Ap,i = p i+ i — p^. Let 



where A N = j^rj. The next lemma shows that the function jF/v is well defined: 

Lemma 3.1. Given N, i, T, pi, AT and Ap it there is a unique set of boundary conditions (Tl, T r , pL 
for the N -chain that leads to these values at site i provided N is sufficiently large (or AT and Api 
are sufficiently small). 

Proof. First, p { and Ap,i determine p L and p R by linearity of p (see Lemma fTTTl ^1.31) , provided Ap^ 
is small enough that both of these numbers are nonnegative. For definiteness, assume AT > 0. 
We will find suitable T L and T R by varying both until correct values are attained at site i: Fix the 
left bath temperature temporarily at T L G (0, T), and vary the right bath temperature from T to 
oo. Since the temperature at site % increases strictly monotonically with the right bath temperature, 
there exists a unique T R = T R (T L ) > T giving the correct value of T at site %. Now look at the 
chain with bath temperatures Tl and T R . As Tl increases to T, AT decreases to 0, so for AT 
sufficiently small, there is a unique Tl for which T i+ \ — T is equal to AT. □ 



C(x) 
C 2 (x,y) 



C(x;T L ,T R , p L , p R ) = 
C 2 (x,y;T L ,T R ,pL,p R ) 



lim N ■ Cov N (S [xN] , S[ xN ]+i), 
= lim N ■ Cov N [S[ xN] , S [yN] ) , 



T(x\ T,p„, Tl, p'J = C(x; T L , T r , p L , p R ) ■ 



(12) 
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With T and F N as above, and for fixed boundary conditions, the limit in the definition of C{x) 
is equivalent to 

lim F N (x; T h pi, ^ J = T(x\ T,,p m ,Tl, pi) ■ (13) 

N->oo y Zijv iiJV / 

Such a limit involves difficult issues beyond the scope of this paper. For example, while the 
convergence of T{ to T* as N — > oo follows from LTE (Assumption 2), the convergence of to 
Tl cannot be deduced from previous assumptions. 

In addition to Assumptions 1 and 2 in Sect. ll .31 we now introduce two other sets of assumptions 
on which the rest of our analysis relies. 

Assumption 3. The functions C and C 2 are well defined and finite-valued. 

Assumption 4. (i) In the context of HI]), ^ — > Tl as N — > oo. 

(ii) The function T defined by 4721) is differentiable. 

Assumption 3 asserts the presence of a well-defined structure at "N = oo" that governs the covari- 
ance relationships of the invariant measures p^ for large N. Such a structure goes beyond the idea 
of LTE to treat information of the next order, namely how p^ deviates from local equilibrium at 
microscopic length scales, and how it deviates from products of Gibbs measures (at varying local 
temperatures) globally. Assumption 4 identifies some technical issues which we take for granted. 

4 Covariances at Microscopic Distances: Pl = Pr 

This section and the next concern nearest-neighbor covariances in long chains, i.e. Cov N (Si, S i+ i) 
for large N. We treat in this section the simpler case where there is no tracer flux in the system; 
the equality pi = p R = p is assumed throughout. Our aim is to discover how the functions C, 
equivalently J 7 , depend on the various quantities. 

4.1 The two middle sites 

In this subsection, we fix x = |, and study the function (T, p, T') i— > F(h T, p, T', 0) where T is 
as defined in Sect. [3] Here T, p and T' are quantities associated with the midpoint of the chain; in 
particular, T' is the temperature gradient at x — |. We begin by examining the dependence of T 
on T' . Fix T and p, and consider 

F(T') =F T >?(T') =F(~;T,p,T',0) ■ (14) 

That is to say, F(T') w N ■ Cgv n (S[n/2], <%v/2]+i) for large N with boundary conditions T L = 
T - \V, T R = T + \T and p L = p R = p. Clearly, F(0) = 0, and F is an even function due 
to the left-right symmetry. This leads one to expect a quadratic response when the system is taken 
out of equilibrium. The following simulation confirms that the coefficient of the quadratic term is 
indeed nonzero. 
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T'(x=l/2) 

Figure 4: Scaling of middle-pair covariance with temperature gradient. The plot shows 
CovjvOSjjV/2]) "SW/^+i) as a function of T" on a logarithmic scale. Boundary conditions are chosen to 
fix T = 60 and k = 2.5^/7r at x = \, Solid discs: N = 8; open circles: N = 20. A solid line of slope 2 is 
shown for reference. See Appendix A for a discussion of simulation parameters. 

Simulation 4. For various pairs ofT and p, compute Govn(S[n/2], S[n/2]+i) far a sample ofT'. 

A set of results are plotted in Fig. H The data show that (i) §fi(0) ^ 0, and in fact, (ii) F(T') 
is fairly well approximated by a quadratic function over the entire range (— 2T, 2T) of T'. (The 
minimum temperature tends to as \T'\ — > 2T.) Note that if F(T') is smooth, then the only other 
possibility is F{T') ~ T' 2k for some k G {2, 3, • ■ ■ }, which is clearly not the case here. Other data 
sets (not shown) with different values of (T, p) confirm these conclusions. 

We investigate next the dependence of a(T, p) : = ^|f&(0) on T and p. First, we establish a 
simple result (which is expected from dimensional analysis): 

Lemma 4.1. Fix T and p, and view A > as a parameter. Define F T,P and F x2T ' Xp as in H4\) . 
Then 

F x2T ' xp {\ 2 T') = A 4 F T '"(T'). (15) 

Proof. Given T, p and T', let T L , T R and p L = pnbe the boundary conditions that give rise to 
these values at the two middle sites, and consider a second system with boundary conditions X 2 T L , 
\ 2 Tr and XpL = ^Pr and initial tank energies A 2 times those of the first system. Via a standard 
coupling argument, the sample paths for these two systems are easily matched. Corresponding 
sample paths give rise to time evolutions with identical tracer counts in both systems - but with the 
second system at energies A 2 times that of the first and running at speeds A times that of the first. 
This leads to ([B]). □ 

Let F x2T ' xp (T') w a x T' 2 . It follows from Lemma|4j]that 

a x ■ (A 2 T') 2 w F A2T ' A "(A 2 T') w A 4 F t ^(T') w AVT /2 . (16) 

Thus a\ is independent of A, and F T,P (T') can be written as F K (T') where k = 2p^pnfT is the 
mean tracer count with respect to p T,p (see Sect. II .21) . The preceding argument tells us that for any 
T, p, k related as above, the function A(k) := a(T, p) is well-defined, and 

Cov N (S N/2 , S N/2+ i) « ^A(k)T' 2 . (17) 
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Figure 5: The A-curve, computed using 8-chains. The boundary conditions are Tl = 20, Tr = 30, and 
PL = PR = 2.5k/ \pn, and A(k) is estimated using Eq. (fTTT ) with N = 8. Left: Linear plot. Also shown 
are some points computed using 16-chains (open circles; Tl = 50, Tr = 100, and pi = PR = 5k/V2tt). 
Right: Log-log plot shows that A(k) ~ 1/k for re » 1. 

A natural question is how A(k) depends on re. Other things being equal, one would think 
intuitively that local thermal fluctuations are smaller when more tracers are present, because stored 
energy is affected less by the entrance of a tracer with very large or very small energy. To confirm 
this, and to collect data for later use, we compute the function A. Practical reasons dictate that 
very short chains be used due to the large number of data points (each one of which requiring a 
separate run) and the time it takes for co variances to converge in long chains; see Appendix A for 
a discussion of relevant numerical issues. However, some understanding of the errors introduced 
by the use of very short chains is necessary if these data are to be useful in later predictions. 

Simulation 5. Compute A(k) systematically for a range of re using chains of 8 cells. Compute 
also some values of A(k) using 16-chains for comparison. 

The 8-chain results are shown in Fig. [5l Note that the function decreases monotonically as ex- 
pected, and A(k) ~ 1/re as re — > oo. 

Let A & (k) and A 16 (k) denote the values of A(k) computed using 8 cells and 16 cells respec- 
tively. We computed A w (k) for some values of re, and find these values to be somewhat larger than 
the corresponding values for A 8 . This is consistent with Fig.[3j and the true A- values (defined at 
"N = oo") are likely to be larger still. We find, however, that the ratios of the two sets of values 
remain fairly constant as re varies. For example, 

1.14 < A 16 (re)/A 8 (re) < 1.20 for k G [2y/n, lOy/n] . 

We will use the 8-chain A-curve data in our study of long-chain covariances. All of our long-chain 
simulations involve k in subintervals of the range shown above, and we will only assume that 
A(k) w const • A s (k). 

4.2 Quadratic responses to local temperature gradients 

We demonstrated in Sect. 14. II that at x = \, the response to the local temperature gradient X" is 
quadratic. We now examine the situation at x ^ \. Recall that our argument for showing that 
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T 

Figure 6: Quadratic response to T' at x ^ \. We fix x / 1/2 and set T(x) = 60, k(x) = 5t/tv. We then 
compute Cova^SVj.tvIi S[xN]+i) as a function of T'(x). Solid discs: x « 0.27; open circles: x « 0.65. 
Here, we use N = 12. We do not take T" much smaller than the values shown for reasons discussed in 
Appendix A. 

the leading term has order > 2 at x = \ uses the left-right symmetry of the chain, a property not 
present at x ^ \. Reasoning physically, however, it is hard to imagine that the response to local 
temperature gradient is sometimes quadratic, sometimes linear, or that the coefficient of the linear 
term would vanish at exactly x — | independent of boundary conditions. The following simulation 
confirms this thinking. 

Simulation 6. For fixed x £ (0, 1), T = T(x) and p, compute Covn(S[ x n], S[ x n]+i) far various 
values ofT' = T'(x). Repeat for various x ^ \. 

Fig. [6] shows the results for some values of x. The function T' i— > T{x\ T, p, T 1 , 0) for fixed x, T, 
and p is clearly quadratic to leading order and contains a nonzero third-order term (which is not 
present when x = |). Other sets of data, computed using different values of T and p, lead to the 
same conclusion. 

Notice that a quadratic response to T 1 is consistent with the observation in Sect. 12.11 that 

Ei^ Cov N (Si, Sj) cx (T R - T L ) 2 ■ N. 

4.3 Boundary effects 

In Sects. 4.1 and 4.2, we showed that (i) for each fixed x, T, and p, T{%\ T, p, T', 0) ~ T' 2 (at least 
for \T'\ not too large); and (ii) at x — |, the coefficient in front of T' 2 is a function of k{\), that 
is to say, T, p, T', 0) ps • T'(i) 2 , the function A(k) being defined by this relation at 

x = \. With A(-) so defined, we now ask if the same relation holds at every x £ (0, 1), i.e., if it is 
true that T{x\ T, p, T', 0) w A(k(x)) ■ T'(x) 2 . 

Consider, for the moment, a different situation in which the system is infinite in length, i.e., 
x E (— oo, oo), and a temperature gradient is maintained by, say, a constant external field. One 
would expect such a system to be translation-invariant, which implies that Cov n (S[ x n], S^tvj+i) = 
Covjv(<SWb S[ y N}+i) if the quantities A(k)T' 2 at x and y are identical, i.e., we would expect the 
answer to the question above to be affirmative. 
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Figure 7: Covariances with identical local conditions at various x. For each x, we choose boundary con- 
ditions so that T{x) = 60, T'(x) = 60, p « 7.75 (closed discs) and 19.4 (open discs). With these local 
conditions, we compute CovAr(SYj,jv]> SLjvi+i) f° r various x symmetrically placed about x = |. Note that 
each point requires a separate simulation. Here, N = 40. We include 1-standard deviation error bars; see 
Appendix A for factors affecting the numerical accuracy of computations. 



Our system has no translation invariance to speak of and, at least for finite N, the randomizing 
effects of the baths are apparent: nearest-neighbor covariances are significantly smaller near the 
boundaries because the baths inject tracers with i.i.d. sequences of energies into the system. The 
answer to the question posed above depends on how fast these randomizing effects dissipate. Do 
they decrease by a certain amount per lattice site, or do they decrease with macroscopic distance? 
The aim of the following simulation is to shed light on these questions. 

Simulation 7. Fix a triplet (T*, p*, X^') and compute Cov n(S[ x n\, S^tvj+i) for various x, using 
(T*, p*, Tl) as the local conditions at x. That is, for each x, we first compute boundary conditions 
which give T(x) = T*, T'(x) = T[, and p{x) = p*, and then use these boundary conditions to 
compute the pair covariance at x. Repeat for various choices of(T*, p*, T^). 

The results, shown in Fig. [71 show clearly that the randomizing effects of baths dissipate in a 
manner more like a diffusion process than by a definite amount per lattice site. Note also the 
asymmetry of the data points: identical local conditions at x and 1 — x do not produce equal 
covariance relations. This suggests that from the conduction point of view, effective distance from 
the boundary may not equal physical distance. 

To capture the phenomenon observed, we introduce, for given boundary conditions Tl, Tr and 
Pl = Pr = p, a function if bc (x) defined by the expression 

C{x) = <p bc (x) ■ A(k(x))T\x) 2 , xG(0,l). (18) 

We add the superscript "be" to stress the fact that this function depends on boundary conditions. 
By definition, <p bc (\) — > 1 as T'(|) — > 0. Since we do not know of any other forces acting on the 
system, we think of <p 6c (x) as a measure of the randomizing effects of the baths at location x, and 
continue our investigation based on this thinking. 
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Figure 8: Comparison of ip bc and ip oc . The boundary conditions here are Tl = 5, Tr = 45, Pl = PR = 8, 
and N = 60. We estimate <p (x) by N • Covn(S[ x n], S[ x n]+i) / As(k)T' 2 where As denotes the 8-chain 
A-curve computed in Simulation \5\ and plot the resulting curve against y = ct(x) where c is chosen to 
minimize the (unweighted) least-squares distance between the two curves. A vertical line is placed at x = \ 
for reference. Clarifications: (i) As noted at the end of Sect. 14.11 there is a number c\ > 1 such that the 
A fa ci • Ag. Thus the data points here represent empirical values of c\(p hc . The constant c\ is absorbed into 
the choice of c, so that the solid curve can be thought of as c\(p hc . (ii) A least-square minimization is used 
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instead of setting the functions equal at x 
visibly unconverged. 



as prescribed in the text because the data near x 



are 



We now seek to identify the function cp (x). To capture the idea of "effective distance to 
boundary," we introduce the following time-to-boundary functions^ in an A^-chain, let t- N ^ denote 
the expected time for a tracer at site i to reach one of the baths. We then have the relation 



T 



r (N) 
i-1 



'i+l 



(19) 



where T, ; is the temperature profile and the expected time a tracer spends at site i. Identifying 

the zth site with Xi = 6 (0, 1) and letting r = -j^t^, we obtain, as A^ — >• oo, the differential 
equation 

t"(x) = tJ== , r(0) = r(l) = . (20) 



Vn 



We define 



<p bc (x) 



where r is the solution of (T20l) and the constant is chosen so that <fi bc (\) 
that ip bc is approximately equal to (p hc 
now test this conjecture numerically: 



const t(x) (21) 

<^ 6c (|), and conjecture 
Note that (p bc also depends on boundary conditions. We 



2 Implicit here is the idea of a dual particle system similar to that for e.g. the KMP model. We do not know if such 
a dual system can be constructed for the random halves model. 
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Simulation 8. Fix boundary conditions and N. For x G (0, 1), estimate ip bc (x) by computing 
Covn(S[ x n], S[ x n]+i) and dividing the result by ^A(k)T' 2 . Compare the resulting function to 

ip bc {x). 

The results are displayed in Fig. [8l They show that to the degree that we are able to estimate 
these functions accurately, the two graphs are in near-perfect agreement. These graphs are not far 
from a perfect parabola, but the definitive presence of asymmetry - in both Fig. 7 and Fig. 8 - 
provides convincing evidence that at comparable distances, it takes longer to reach the bath at the 
lower-temperature end of the chain. 

Remark. When the system is in equilibrium, i.e., when Tl = Tr = T, the solution of (T20l) 
is easily computed to be t(x) = 775^(1 ~ x )- Recall that the expression x(l — x) also appears 
in the formulas for spatial covariances in the simple exclusion Q|20l and KMP models [3]. An 
important difference between these models and ours is that there, the clocks signaling state changes 
ring at rate 1, whereas in ours, they ring at energy-dependent rates. This property is responsible 
for, among other things, the slight asymmetry in r. 

From the limited data available, we clearly cannot conclude the validity of our conjecture, but 
details aside, the general ideas seem to point in the right direction. 

5 Covariances at Microscopic Distances: Pl^ Pr 

We return to the general case where both energy and tracer fluxes may be present. The situation 
here is more complex, and we are less able to separate the contributions of the various factors. We 
will focus on the dominant effects, and make some observations along the way. 

5.1 Responses to local gradients 

As before, we consider first x =\. For fixed T and p, we consider the function 

F(T', p') = F T >»(T', p') = T, p, T\ p') 

and begin with a few observations on which terms may be absent in its Taylor series. 
Proposition 5.1. F(0, p') = Ofor all p'. 

Proof. Using the notation in Sect. [T] we claim that when T' = 0, the invariant measure Pn of the 
chain is the product 

p T ' pl x p T ' p2 X • • • x p T ' pN 

where Pl = P — \p' and Pi = Pl + jn^p' for z = 1, • • ■ , A^, so that stored energies at distinct sites 
are uncorrelated. This is proved by direct verification (along the lines of the proof of Proposition 
3.7 in |[9]|). See Appendix B for details. □ 
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Figure 9: Local polynomial expansion of F. We set T = 100, p G {12.5, 20}, T' = t, and p' = r at x = ^, 
then compute F(t, r) ~ N ■ Covat(^) for N = 8 using the (t, r) values from the grid (left panel). We then 
find the best polynomial (in the sense of least-squares) of the form at 2 + btr + ct 4 + dt 3 r + et 2 r 2 + ftr 3 
which fits the data (right panel). For both values of p, the absolute standard error for the coefficient of the 
tr term is ~ 2 x 10 -3 . 



Interchanging Tl with Tr and pl with pr, we see that F(T', p') = F(—T', —p'), so that there 
are no odd order terms in the polynomial approximation of F(T', p'). Proposition 15 . 1 1 says there 
are no terms that are purely powers of p', and we have shown in Sect. |4.1| that among the terms that 
are pure powers of T' , the leading one is T' 2 . 

Perturbing from equilibrium, we have shown that ^-(0,0) = |^(0,0) = 0. Consider the 

Hessian D 2 F(0, 0). Since = (Proposition O), it follows that for \T'\, \p'\ < 1, 

F(T', p') w aT' 2 + bT'p', for some a = a(T, p) and b = b(T, p). (22) 

Recall our brief discussion of the sign of Covtv(-, •) in Sect. 12.21 Our next lemma says that the 
degeneracy of the Hessian D 2 F(0, 0) is equivalent to all covariances at x — \ having the same 
sign, namely the sign of a(T, p). 

Lemma 5.1. F > in a neighborhood of (0, 0) if and only if a > and b = 0. 

Proof. As a quadratic form, the Hessian D 2 F(0, 0) has maxima and minima in perpendicular 
directions. We know already that it is degenerate along the //-axis. Thus either the minimum is 
negative, that is F(T', p') < for some (T", p'), or the Hessian is positive semi-definite, i.e. a > 
and 6 = 0. □ 

To determine the leading-order terms of the Taylor expansion of F near (0, 0), in particular to 
see if b = (equivalently F > 0), we perform the following simulation: 

Simulation 9. Fix T and p, and set T' = t, p' = r. Compute F(t, r) for a grid of (t, r) values and 
use the resulting data to find the best fitting 4th degree polynomial of the form 

Pit, r) = at 2 + btr + ct 4 + dt 3 r + et 2 r 2 + ftr 3 . (23) 
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Figure 10: Responses to local gradients for various T' and pf at x = |. We fix T, p, and T' at x = ^ 
then plot N ■ Cow n {\)/A{k)T' 2 as a function of p'/P- Each curve corresponds to one value of T 1 . The 
parameters are N = 12, T(±) = 50, p(~) « 9.68, and T'(±) = 90, 70, 50, and 30 (top to bottom). Some of 
the curves are discontinued because when T'(^) and p'{\) are large, T(x), which decreases monotonically 
with x, may reach before x does. This is easily seen from TheoremQ] 



Two sets of results is shown in Fig. [9J We find that in these simulations (and in others not shown 
here), the leading-order term in the Taylor polynomial of F(T', p') is aT' 2 with a > 0. In all cases, 
\b\ <C 1, which is consistent with all observed covariances being > 0. As expected, we cannot 
conclude definitively that 6 = 0. 

From here on we take as a working assumption that the quadratic form associated with D 2 F(0, 0) 
above is degenerate, i.e., at x = |, the only second order term is A(k)T' 2 as in the p L = p R case. 
For similar reasons as before, namely that the qualitative picture elsewhere in the chain should not 
differ from that at x — |, we assume further the absence of T'p'-terms in the Taylor expansion of 
F at all x E (0, 1). This completes our discussion of how covariance depends on local quantities 
when the system is not far from equilibrium. 



5.2 Far-from-equilibrium corrections 

We have shown that at x — \, F T ' P (T', p') w A(k)T' 2 is a good approximation of the response 
to the local temperature and injection gradients - provided \T'\ and \p'\ are sufficiently small. In 
the p' — case, we noted in Sect. 14. II that the function T' i— > F T,P (T') is close but not equal to a 
perfect quadratic for large T' . We now investigate the effect of larger \p'\ on F T ' P (T', p'). 

Simulation 10. For fixed values of (T, p), compute F T ' P (T', p') at x = \as a function ofT' and 
p' , and study the deviations from the function A(n)T' 2 . 

A set of results is shown in Fig. |T0l The results for other choices of (T, p) (not shown) are similar. 

Observe that (i) the dominant factor is A(k)T' 2 , but a moderate correction (roughly 20% for 
\p'/p\ < 1) t0 me leading coefficient is sometimes needed; and (ii) a larger \p'\ tends to increase 
covariances. Along the lines of the conjectural thinking that positive covariances are caused by 
waves of abnormally high (or low) energy tracers with long wavelengths, the presence of a tracer 
flux appears to amplify the propagation of these waves. 
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5.3 Nearest-neighbor covariances in long chains 

We have seen that the presence of a tracer flux complicates the behavior of covariances at x = 
\ when the system is far from equlibrium. We do not know how it affects the propagation of 
boundary effects, or how to separate these contributions (or if they can be separated at all). We 
will demonstrate in this subsection, however, that the main ingredients in Cov n(Su Si+i) are those 
already identified in Sect. |4] 

More precisely, let (p bc and (p bc be as defined in Sect. 14.31 That is to say, ip bc is as defined in 
(TT8T) . r the solution of (1201) . and (p bc (x) = ct(x) where c is chosen to ensure cr(\) = <£> 6c (|). We 
test the validity of 

Cov N (S [xN] ,S [xN]+1 ) « 1 • bc (x) ■ A(k(x)) ■ T'{x) 2 (24) 
as an approximate relation. 

Simulation 11. For fixed boundary conditions with pi ^ p^, and N taken as large as possible, 
compare empirically computed values of Cov tf(Si, S i+ i) to their predicted values given by d2?1) . 

The results for two sets of boundary conditions are plotted in Fig. \TT\ Notice that only the shapes 
of C(x) are being predicted because of unknown normalization factors. 

Fig. [TT1 confirms many features of our predictions. First, the shapes of C(x) confirm that the 
randomizing effects of the heat baths are diffusive in nature; these effects result in an effective- 
distance factor whose graph takes the form of a distorted parabola. Second, the locations of the 
peaks of C(x) show the influence of T'\ see the left column of Fig. 1 1 HI Finally, comparison of the 
data for different N in each set provide another confirmation of the 0(1/ N) scaling. 

Fig. [TT] also shows the need for corrective factors on the order of 10 — 15% away from the 
boundaries. While we do not know the precise nature of these corrections, we point out that in 
the presence of a tracer flux, the functions tp bc (as computed from empirical data) are a little more 
asymmetric than in the case pi = Pr (shown in Fig. [8]). 

6 Long-range Covariances 

6.1 Renormalizability and the function C2(x, y) 

In Sect.[3l we introduced the idea of a pair-covariance function C 2 (x, y). The existence of the limit 
in the definition of C 2 (x, y) implies the following: for all sufficiently large N, N', 

N ■ Cov N (S[ xN ], S [yN] ) w N' ■ Cov N '(S{ xNI ],S{ yN/ ]) . (25) 

This can be seen as a statement about the renormalizability of pair covariances: Consider an iV- 
chain with N = rN for some integers iVo and r, and subdivide the chain into iVo groups of r 

3 In Fig. 11(a), both T' 2 and A increase from left to right; in (b) T' 2 increases from 4 x 10 2 to 6 x 10 3 while A 
decreases from 0.07 to 0.04. 
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(a) T L = 10, T R = 50, p L = 15, p R = 9; N = 70, 90, 120 




(b) T L = 5, T R = 45, p L = 8, p R = 15; N = 50, 70 



Figure 11: Nearest-neighbor covariances. Left column: Temperature profiles T{x). Right column: Em- 
pirically computed values of CovjvOSi, Sj+i), plotted against predicted covariances using Eq. (|24l) . The 
predicted covariance curves are scaled to coincide with the empirical curves at x = \. 

consecutive sites. For convenience, we take r to be odd. Let = (i — 1)N + |(r + 1), i.e., 
is the index of the middle of the r sites in the ith group. For any i, j with 1 < i, j < N and i ^ j, 
we compare the covariance at sites i and j in an iVo-chain to that at sites and j( r ' in an iV-chain 
with the same boundary conditions. Eq. (|25T) tells us that r ■ Cov r jv (^M , S^-w) ~ CoVjv (<Si, Sj) 
if N is sufficiently large, and that as r — > oo, r ■ Cov^^m, converges to a constant. 
We test this renormalizability to confirm that C 2 (x, y) is well defined. 

Simulation 12. Fix T L , Tr, Pl, Pr, arcd iVo G Z + . For various pairs in the N -chain, 
compute N ■ Cov A r(S , i (r) , Sj( r ))for N = rN , r = 3, 5, 7, • • ■ . 
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Figure 12: Renormalizability of the covariance function. The plot shows N • Covj^(S^r) , S^)), where N = 
tNq, Nq = 8, and r = 1, 3, 5, • • • ,13. From top to bottom, the curves correspond to the (i, j) pairs (3, 4), 
(3, 5), (2, 3), (2, 4), (2, 5), and (2, 6). Boundary conditions are T L = 10, T R = 100, p L = 20, p R = 10. 




Figure 13: Covariances at macroscopic distances. These plots show Covn(S^ xN ^, S^ yN ]) as a function of y, 
with x fixed at ps 1/3. 

A subset of the data is shown in Fig. [121 Because the value of iVo used is relatively small (No = 8), 
one can expect the plotted values rN ■ Cov N (S i ( r ), S^-m) to converge or stabilize to constants only 
as r increases. The results show that they, in fact, stabilize fairly quickly. 

That C 2 (x, y) is well defined implies that the function (x,y) i— > Cov N (S[ xN ], S[ yN ]), when 
normalized, settles down to a fixed shape for large N. We now investigate the shapes of these 
functions. While carrying out Simulations [8] and [121 we also collected data for Cov^(Si, Sj) for 
various pairs of Graphs of j ; i— > N ■ Cov^(Si, Sj) with i = [|iV] are shown in Fig. [13] for 
three sets of boundary conditions. As predicted in Sect. 12.21 these functions are bounded, and they 
decrease monotonically to as \i — j\ increases. This decay rate is roughly linear in macroscopic 
distance, and extraordinarily slow per lattice site. For example, if N = 10 6 and i = [|A r ], then 
Cov^iSi, Si+i) ~ 2 • CovAr(5'i, S21). We note that these findings are consistent with fluctuating 
hydrodynamics; see e.g., [|20ll . 
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Figure 14: A slice of the 2-point distance function ip\f. Here, we use Eqs. (l26l ) and (|27T ) to compute the 
graphs of y i— > ^(x, y),x « |, from the data in Fig. [13] The results are normalized to 1 at x and overlaid. 



6.2 An approximate formula 

A natural generalization of the nearest-neighbor covariance formula (1241) to pair covariances sepa- 
rated by macroscopic distances is 

Gov N (S [xN] ,S [yN] ) « — -(pf(x,y) ' A(x,y) -T'(x,y) 2 . (26) 

We think of this as an approximate formula that holds for x,y G (0, 1) with < \x — y\ <C 
1: ip b f(x,y) is a notion of effective distance of the pair x,y to the boundary, and A and T" are 
generalizations of corresponding quantities in d24l) . As a rough approximation, one may take 



A(x, y) = \ (A(k(x)) + A( K (y))) , f '(*, y) = T(y) T(x) . (27) 
2 v 7 y — x 

To obtain information on (f\ c (x, y), we estimate it using (1261) . taking A and T' as above. Slices 
of the graphs of these functions of 2 variables, with x fixed and y varying, are shown in Fig. IT4l 

Piecing together these slices, we deduce the following geometric facts about our 2-point dis- 
tance function Lp b 2 c defined for (x, y) E (0, 1) x (0, 1): it is continuous and piecewise smooth, with 
a "ridge" along the line {x = y}. Along that line, it is equal to p bc (x) = (p hc (y), the value of 
the 1-point distance function introduced in Sect. 14.31 For fixed x, the function y i— > Lp\ c (x,y) is 
roughly piecewise linear, peaking at y — x. 

Observe that if this function were exactly piecewise linear, then by exchanging x and y, we 
would arrive at the relation 

V bc (x) ■ = v" c (y) • - for x < y . 

1 — x y 

This relation implies that ip bc (x) = cx(l — x) for some constant c. Since (p bc is roughly parabolic 
(see Sects. 1431 and 15.31) . the approximate piecewise linearity of <^| c is consistent with our results 
from previous sections. (Exact formulas for the simple exclusion [01201 and KMP models [3] also 
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contain piecewise linear functions.) On the other hand, we know from Sections 4 and 5 that for the 
random halves model, ip hc (x) ^ cx(l — x). Thus ^ cannot be exactly piecewise linear either. 

In summary, our results for long-range covariances are rough and are obtained by extrapolat- 
ing from what we know about covariances at microscopic distances, together with a numerical 
determination of the distance function An approximate formula is 

i f -(i5f)-^(x)-A(x,y)(^ M ) 2 for 0<x<y<l, 
^.(j£).^(x).ii(x,y)(S^5i) a for < y < x < 1 , 

where (p hc is our approximate one-point distance function and A is as in (ITTT) . Both this formula 
and Fig. [13] show clearly the following: 

(i) Covariance decays essentially linearly with macroscopic distance, i.e., 

C(x) -C 2 (x,y) ~ \x-y\ . 



(ii) For fixed x, the function y i— > £2(2, y) is continuous but not differentiable at y = x; it has 
opposite concavity on the two sides of x if the temperature profile T(y) is nonlinear. 

That the curvature changes sign at x is clearly visible in Fig. [T3l it is also evident from the second 
derivative of (IM^He1\ 2 ( see Fig.QTl left column, for the temperature profiles T(y)). 



Conclusions and Remarks 

Via a series of theoretical arguments and numerical simulations, we have developed a coherent pic- 
ture for the spatial covariances at steady state of the 1-D random halves model. We have established 
firmly that stored-energy covariances have order of magnitude 4 away from the boundaries of N- 
chains. This in itself points to the presence of long-range covariances which decay very slowly. 
Subsequent analysis shows that covariances decay linearly with macroscopic distance. For sites 
separated by microscopic distances, we have a simple formula that encapsulates the main ingredi- 
ents on which energy covariances depend, including (i) a quadratic response to local temperature 
gradient, (ii) diffusive nature of the randomizing effect of the heat baths, and (iii) stabilizing effects 
of large numbers of tracers. 

Since the random halves models are stochastic idealizations of certain mechanical models, we 
hope our results will also shed light on these and similar Hamiltonian systems. There are, however, 
important differences, such as mixing and memory issues (see e.g. [|9l[8l). The extent to which the 
picture established here carries over to the Hamiltonian setting remains to be seen. 

Finally, it is well known that nonequilibrium phenomena are quite different in higher dimen- 
sions. Our detailed study here provides a baseline for explorations in 2-D and 3-D, a project 
currently being carried out by the authors. The results, which are indeed different as predicted by 
fluctuating hydrodynamics, will be reported in a forthcoming paper. 
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Appendices 

A Remarks on simulations 

The simulations carried out in this paper implement directly the dynamics described in Sect. Il.lt 
expectation values with respect to the invariant measures are computed via time averages over 
long trajectories. The numerical issues are quite similar to those of Markov chain Monte Carlo 
computations. 

We calculate empirical error bars using a standard "batch means" estimator. These error bars 
measure only statistical errors that arise from the fact that expectation values are estimated by time- 
averaging over finite intervals of time, i.e. finite-length trajectories. The error bars do not account 
for finite-size effects, i.e. bias due to the fact that we can only simulate iV-chains with finite N. To 
improve clarity and readability, we have suppressed the error bars in most of the figures, displaying 
only those that are directly relevant to the issues being discussed. 

A variety of variance reduction techniques have been invented to speed up the convergence 
of Monte Carlo calculations, ranging from multigrid Monte Carlo |[T9l to large deviations-based 
importance sampling [4J. Most of these techniques require additional information, such as an 
explicit expression for the invariant measure or a large deviations functional. Such methods do not 
work in our setting. There is on-going work on a class of algorithms which do not require detailed 
knowledge of the invariant measure |[T2|. However, it is not known whether such techniques can 
be applied to the random halves model. 

The quantities of interest in most of our simulations are the Cov N (Si, Sj). In general, these 
are relatively small numbers that are differences of two much larger numbers and can be rather 
costly to compute. We estimate the cost of computing Covjv as follows: for temperatures which 
are ~ T, both E(S i S j ) and E(5' i )E(S i ) are ~ T 2 , while Cov N (S h Sj) ~ A(k)T' 2 /N, where R 
is the typical tracer density per site. (For the present discussion, we focus on i,j away from the 
boundaries so that (p bc (x) ~ 1.) In order to compute Cov N (Si, Sj) with a relative error of e, we 
need E(SiSj) with a relative error of eA(k)(T' /T) 2 / N; the same is true for E(S i )V l (S j ). Since 
the statistical error in the time average i J Q T Si(r')Sj(r') dr' is ~ \fotfr for some a, this means 
we need to integrate the system for a time r which is proportional to N 2 /e 2 A(R) 2 ■ (T/T') 4 . Now, 
the computational cost of simulating the system up to time r, as measured by the total number of 
"events," is proportional to the number of tracers RN in the system and the mean rate of activity of 
each tracer (which is ~ VT). Thus, we have 



For our simulations involving long chains, we have found a to be typically < 10. As an example, 
for e = 5%, N = 60, T L = 5, T R = 45, p L = p R = 8, we have R w 6, and A(R) ~ 0.1, so that 




(28) 
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~4x 10 11 events are needed. On our computer system0this requires ~ 14 days. 

Eq. (|28l ) tells us that the computational cost grows rapidly as T'/T decreases. In Simulations^] 
[6l and [91 where we study the small-T' behavior of the covariance function, this rapid growth 
prevents us from taking T' too small. However, we do not always need small T'/T, and a large 
temperature gradient not only reduces the computational cost, it also takes the system farther out 
of equilibrium so that some effects are made more transparent. 

The cost also grows rapidly as N increases. This is why we use short chains wherever possible, 
for example in computing the A-curve. In Simulation [5j we use N = 8 because we feel that the 
amount of accuracy gained from using longer chains is perhaps not worth the additional cost. 
However, it is not always possible or advisable to use short chains, particularly where infinite- 
volume limits are involved. In such cases, one can sometimes obtain better results by increasing 
T'/T. 

Finally, recall that A(k) ~ 1/R for R ^> 1, so the computational cost is ~ R 3 for large R. 



B Invariant measures 

Proof of Proposition \5. 1 [ The proof of invariance of the measure fi T ' pl x • • • x n T < PN follows closely 
that of Proposition 3.7 in [9|. We refer the reader to flU for some of the background notation. Fix 
a phase point 

z = (z« = ({x?,...,^},v il] ; {x[ N \...,xl% y W). 

We assume x± , . . . , x^ are distinct, and let e, h > be arbitrarily small numbers. Fix arbitrary 
n with 1 < n < N . We will compare the 3 probabilities, P n ., P n+ i >n and P n _i jn defined below: 

Let P n> . be the probability that at time t — 0, in every site j, there are exactly kj tracers the 
energies of which lie in disjoint intervals 

[x?,x? + e], [xg>,xg>+ e ], 

and before t = h, the tracer in site n with energy in [x^\x^ + e] exits the site. The number 
P n +\,n is the probability that at time 0, the tracer configuration is as above except that the tracer 
with energy in [x\ , x^ + s] is in site n + 1 instead of site n {i.e., there are k n — l tracers in site n 
and k n+ i + 1 tracers in site n + 1), and before time t = h, this tracer jumps from site n + 1 to site 
n. The number P n -\,n is defined analogously with site n — 1 playing the role of site n + 1. 

To prove the invariance of /i T ' pi x • • • x /j T ' PN , three sets of balancing conditions have to be 
met. The interaction with tanks is as in [EY]. We verify below the equation involving interaction 
between neighbors, namely P n> . = P n +i : n + P n -i,n, an d leave the one involving interaction with a 
bath to the reader. 

4 We performed most of our simulations using GCC on 900 MHz SPARCv9 processors. Note that this running time 
estimate also depends on how much information is collected during the simulation. 
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Let Ofc, cy and py be as in Proposition |l.l| (Sect. |1.2| ). We will use the shorthand a ki = a ki ( z ^)^ 

(i) 

and write p k to remind ourselves that p = pi at site i (c^ and afc do not depend on p). First, 



^(n) 



(n) 



where Z is defined by the equality above. Next, -P n +i in = \l ■ II ■ III where 

/// = p£tJiiC iWl+ ia iWl e*- +1+1 f* "^^e-^x. 
This product can be written as 

h z f Pk^-i c kn -i \ / Pfcnti+i c fcn+i+A 



2 v p£ } c *» J v pK } / 2 V2p„ 

The equality above follows from the relation 

CkPk T_ 
c k+1 p k+1 2p ' 

which can be derived from the characterization of p T ' p (Proposition ll.il) . An analogous argument 
holds for P n -i, n , an d the desired equality follows from 

1 [ Pn+l , Pn-\ \ ^ 



2 V Pn Pr. 

□ 
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